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Abstract. We consider a smooth, rotationally invariant, centered gaussian process in the plane, 
with arbitrary correlation matrix C^^i. We study the winding angle <j>t around its center. Wc 
obtain a closed formula for the variance of the winding angle as a function of the matrix C^^/ . 
For most stationary processes Cui = C'(t — t') the winding angle exhibits diffusion at large time with 
diffusion coefficient D = dsC" (s)'^ / (C'{0)^ — C(s)^). Correlations of cxp(in(j>t) with integer n, 
the distribution of the angular velocity tpt, and the variance of the algebraic area arc also obtained. 
For smooth processes with stationary increments (random walks) the variance of the winding angle 
grows as ^(Int)^, with proper generalizations to the various classes of frax;tional Brownian motion. 
These results are tested numerically. Non integer n is studied numerically. 
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1. Introduction and model 

The winding of planar random processes has been studied for a while. These are of interest for the 
physics of polymers [TJ [5] , flux lines in superconductors |3] H] and quantum Hall effect [S] . Recently 
there was revived interest in winding properties of processes described by Schramm-Loewner Evolutions 
(SLE) [6], such as the loop erased random walk j7]. In the case of the planar Brownian motion the 
distribution of the winding angle 0t around a point O was computed a long time ago by Spitzer [8] 
who found that t/fj/lnt has a Cauchy distribution, i.e with infinite first moment. This peculiar feature 
was later understood to be related to the large winding accumulated while the trajectory wanders 
infinitely close to point O, and is removed by considering either a small excluded region around O, or 
a lattice cutoff, or some other regularization of the short time, and leading instead to exponentially 
decaying distributions [9j [lOj EJ [12] . Recently it was found that correlated gaussian processes related 
to the fractional Brownian motion, a scale invariant process with stationary increments correlated in 
time, has similar properties [IB] . 

The aim of this paper is to study the winding of a very general continuous-time gaussian process 
= + i^^ in the complex plane with arbitrary correlations in time. The only restriction, mainly to 
avoid cumbersome formula, is that the measure is rotationally invariant around the origin and the 
winding angle is measured around point 0, i.e. = rje*"^' where (jit is a continuous real function 
of time and = \£^t\. The process is thus centered (^j) — and fully characterized by its two-time 
correlation function: 



with i,j = x,y, equivalently {^tS.t') — and (^jCf) = '^Cw- The most general form would be 
(CtCf) = 2(C't(' +iAtt') but we also assume reflection symmetry which forbids the antisymmetric term 
Cy^tt' in (Uj) with ei2 = —£21 = 1. Since the process £,t/\/Ctt has the same winding angle as ^t, 
observables involving only the winding angle should only depend on the combination: 



with Ctt — 1 and from Cauchy-Schwartz inequalities, \ctt'\ < 1, the bound being saturated, i.e. 



Ctt' — ±1, if and only if — ±y ^t'- Some particular cases are (i) stationary process 
Ctt' — C{t — t') and one defines Cto — c{t) — C{t)/C{0) (ii) process with stationary increments 
Cjj','"'^' :— dtdfCtt' — C2(t — t') (here and below wc adopt the following definition for partial derivatives 
dtCtt' — C^l/^^ etc.). (Strict) normalizability of the Gaussian measure requires that these functions 
have (strictly) positive Fourier transforms c{uj) > and C2{oj) > 0. In some cases we restrict to 
a process which everywhere below we call "smooth", meaning - by definition here - £,t differentiable 
at least once, i.e. C^/'^'' exists (equal — C"(0) for a stationary process). For such a smooth process, 
iit^t) — Cfl'^'^ = C'tt''^'' = \dtCtti which vanishes if furthermore the process is stationary. These 
processes appear everywhere in physics, prominent examples occur in the study of quantum noise, 
where C{u)) = r^huj cot\i{hu) /2T) (supplemented by a bath-dependent short time cutoff, e.g. at T = 0, 

C{lu) — firy|w|e~'^°l"l , i.e C(t) = j^S^^ryi) or of quantum Brownian motion [14 and has served as a 
motivation for the present work [TS] 

The outline of the paper is as follows. In Section 2 we study single time quantities. The distribution 
of angular velocity is obtained. In Section 3 we study the periodizcd winding probability distribution 
which is easier than the full one. The correlations of exp(in(j)t) are obtained analytically for integer n, 
and studied numerically also for non-integer n. In Section 4 we obtain a closed formula for the variance 
of the winding angle as a function of the matrix Ctt' ■ We show that for most stationary processes the 
winding angle exhibits diffusion at large time and we obtain the diffusion coefficient. We also study 
non-stationary process such as the random walk and the various classes of fractional Brownian motion. 
Finally in Section 5 the variance of the algebraic area is obtained. Most results are tested numerically. 



i^et') = k.Ctt' 



(1) 




(2) 
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2. Single time quantities 

Single time quantities are easily extracted from the Gaussian distribution ~ d^^^e"'^''^/^^*^"^ 
performing change of variables. Everywhere below (P^ = d^d£^* = d^^d^''. The modulus is 
distributed as P{rt)drt with P{r) = /{"^Ctt) ^ hence the probability to be within rt < t near 

the center vanishes as e^/(2Ctt). To compute the distribution of the angular velocity 4)t one uses 
that Xt — {(,t,it) is gaussian with measure ^^-^^^^p^ det(Af )e^3^ and correlation matrix M^^ = 
{{Ctt,cj:l'^^), (Ct°'"^\ Cjj"^'^^)). Let us denote — at£,t with a* = ft/r* + i4>f Here we have requested 
a smooth process. The measure becomes ^'-^^j#^|^tpdet(M)e~5^l^*l^ where (3 = (1, )M(1, at). 
Integration over yields the joint distribution P{pt, (f>t)dptd(j)t, with pt = Inrj, equal to: 



dptd(pt CttCtt 



(1.1) 

(3) 



Integration yields: 

P(0O4t = 4. , (4) 
2{at + (gy/^ 

with at — {CttCfl'^^ — {Cj:l'^^)'^)/Cff = dtdt'ln\ctt'\\t'=t- For a stationary process at = a = — c"(0). 
For stationary increments at — C2{0)/Ctt — \ {dtlnCtt)'^ ■ Note that this distribution is broad, it does 
have a first moment but no second moment i.e. {(f)t) is infinite. 

3. Periodized winding 

Next one can compute two time correlations of the winding angle. The two time probability measure 
of the process can be written: 

rtrt>drtdrt'd(l)td(f>t' , Cft'rt + C'ttrl ~ 2Ctt'rtrt' cos{(f>t - (l>t') . 
i2n)^At, ) 

with Aft/ = CttCt't' "C't^f hence integration over rt and rt> allows to obtain the probability distribution 
of cos((/)t — (f>t')- Equivalcntly this gives the probability of (p :— (pt — 0t' modulo 27r, i.e it gives the 
periodized probability P(0) = X]m=-oo + 27rTO) where P(0) is the probability of the total winding 
g] — cxd, +cxd[. Defining: 

^, , [ _x2_y2 \ zarccosf— 

F{z) ~ / xye ^ = ^ + 6 

Jx,y>Q 1 - ^ (1 - 

for -1 < z < 1, with arccos(-z) ^ 7r-arccos z and = ^^^^^ |j-r[l + f]2z" = l + 2^ + 2z2 + 0(z3), 
one finds: 

P{cP)^^{l-cl,)F{ctt'COsm (7) 

One can check that J^^ d(j)P{(f)) = 1. An interesting limit is t close to t' . Then ctr is close to unity 
and using the expansion F(z) — „ , ^ + 0(1) we find that: 

and furthermore we expect that ~ ^"(0) since the probability of 27m winding with n 7^ is 

negligible in that limit. Upon Taylor expansion in i — t' one finds that this result is consistent with ^ 
but the approach here is more general as the formula ([t]) does not require a smooth process. The only 
assumption in ([8| is then the continuity of Ctt' ■ The short time behaviour is controled by the distribution 
for moments with < a < 2, i.e. p « K^^il ~ Ctf)"'^^ with X„ = 2"/2r(l - |)r(i±^^)/0F, 



Winding of planar gaussian processes 



4 



and becomes dominated by the cutoff at </> = 0{\) for a > 2 with (|0|")p ~ K'^{1 — Cu'Y ^- For 
instance the variance of the winding angle is found as 

(02)p«-(l-QtOln(l-c«O (9) 

for Ctt' close to unity, and we check below that this coincides with the behavior of {4>'^)p at short time 
differences. 

This allows to compute the correlation functions Cn{t,t') — (e*"'-''^'"'^''') for integer n, which thus 
have closed expressions as a function of the matrix Ctt' '■ 

C^{t,t') ^ F^ictt') (10) 
One finds for instance: 

F^ic) = -^iEic') + {c'-l)Kic')) (11) 

F,ic)^l + {\~l)Hl~c') (12) 
and we obtain the limiting behaviours: 

^^i(c) = ^ + ^ + ^ + 0(c^) (13) 

= 1 - ^(ln(-^) - 1) + 0(ln(l - c)(l - cf) 
1 1 — c 

for small c (large time separation), and for c near 1 (small time separation), respectively, and 

F^ic) = ^ + 0(c"') and more generally F„(c) = -^-^Tp^c" + 0(c"+^) for integer n at small c. 

We have checked these results numerically for several stationary processes where Ctv = c(r) = 
C(t)/C(0) where r = t — t' . The process Q was generated numerically using a discrete Fourier 
transform of -y/c(w)A^Ar^% where N the number of points is typically N = 2^^, At — .01 is the time 
segment in the process and is a unit white gaussian process. We have computed C„(r) where the 
average (e™*^) is over the time range and over several realizations, typically 10. We have plotted Cn(r) 
parametrically as a function of c(r) for various type of noises. Up to numerical accuracy all the curves 
fall on the predicted master curve C„(r) = i^„(c(T)). When c(t) is non monotonous, the master curve 
may be traced more than once. This is illustrated in Fig. [T] 



4. Variance of the total winding angle 

The previous results are easy to derive, and are simple functions of Ctt', but they do not contain 
information about integer winding. They only probe P{4>), the periodized winding angle distribution. 
An interesting question is how to access the full winding distribution P{(f>) and whether its dependence 
on the matrix Ctt' remains tractable. It is a more difficult question since to compute the full winding 
angle one must follow somehow the time evolution of the process, e.g. use that cj) = (pt — 4>t' — J^, dcfig. 
A related difficult question, which requires the full distribution P{4>), is to obtain the averages 
Cn{t,t') = (^e'"^^'^t-<Pt')^ for non integer n. It is seen on Fig. [2]that these are not simple functions, but 
rather unknown and more complicated functionals, of Ctt' ■ 

Here we present the simplest result on this question, the variance of the winding angle. Here, for 
simplicity, and to avoid stochastic calculus subtelties, we first restrict to a smooth, i.e. differentiable 
process as discussed above. We only need to compute the two time angular velocity correlation 
C^{t,t') = {Mt')- We use that ^t = Im(^r)/ICI' = ^^J^lill\^t? hence: 

C,{t,t') = 2[ Gss'tt' (14) 

J s,s'>0 

Gss'tt' — {it^t')s,li'{it^t')s,s' — {£,t^t')s,s'{^tQ)s,s' 
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c(t) 

Figure 1. The correlation function C„(r) for n = 1 (top curve) and for n = 2 (bottom curve) as 
a function of c(t) for c(t) = (Blue curves) and the prediction for -Fi(c) and F2{c) (Red 

curves). The inset shows how the curve is retraced for negative values of c. 




c{r) 

Figure 2. The correlation function Cn(T) for n = ^ (three curves starting on the top) ^ (three 
curves starting on the bottom) as a function of c(t) for three processes: c(t) = exp— (in blue), 
the c(r) used in Fig.l (in green), c(t) = (1 — T'^)exp— (in red). Note that for each n the three 
curve remain very close for c > 0.4 and that for n = 5/4 all processes change sign 
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using isotropy, where {0[^x])s,s' — (0[^2:]e~''^^'^^^~*'^*t''^^) and the same for averages over £_y. The 
integrals over s restore the (difficult) denominators. Each average can be computed from the generating 
function: 

y, g-A'lM2Ctjt2-Z;i = i,2 3C'titiA'i+A'i(A'Ctit+M'Ct.j/) ^^5) 

by taking appropriate derivatives w.r.t. and at /x^ = and coinciding times, and integrating with 

the measure ^^^^ ^, e~^~ to restore the (...)s,s' averages. After a straighforward but lengthy 
algebra, and massive simplifications, one finds: 

^(1,0)^(0,1) , ^(14)^ 



(1 + 2sCu + 2s'Ct.t. + 4ss'(C„Ct-t. - C^,))' 
which, after integration over s, s' gives the angular velocity correlation, which is our main result: 

1 ^(1,0)^(0,1) 

= -^(at9*,ln|c,H)ln(l-4,) (16) 
where we recall that cu' — Cu' I \JCttCt't' ■ The variance of the winding angle is then obtained as: 

$tt' = ((0* - 4>t'f) = rftl^^ dt2Cy{h,t2) (17) 



hence dt^tt' = 2 /j! dt2Cv{t, t2), where Cy{ti, t2) is given by ( 16 1. We now discuss separately stationary 
and non-stationary processes. 

4-.1. stationary processes 

Let study first stationary processes Ctt' — c(r) — C(t)/C(0) with t = t — t'. Then the angular velocity 
correlation becomes Cy{t, t') = Cy{t — t') with: 

C.(.)4(^™^)ln(l-c(.)^) 

^-(a^ln|c(r)|)ln(l-c(r)2) (18) 



2 



which exhibit a divergence at small time r, C„(t) « c"(0) ln(T^-c"(0)), but an integrable one. The 
winding angle variance takes the form ^u' = ^{t — t'), and using that ^^-^(r) = 2 JJ" dT2Cy{T2) one 
finds upon integration by part: 

9,$(r) = 2 r ds-^^^^ + ^ ln(l - c(t)2) (19) 







1-C(s)2 c(t) 

with no boundary term at r = since c(r) w 1 + ^c"{0)t^ at small r, with c"(0) < since 
the process is smooth. The small r behaviour of the variance of t he winding angle is <&(t) 
T2(c"(0)ln(Tv/-c"(0)) - |c"(0)). More generally, the formula (19 1 holds for processes such as 



c(r) = e~'^ at small r with a > 1 so that the small time singularity of C„(r) be integrable. 

If we now consider processes such that c(+oo) = then we find that the generic behavior is that 
the winding angle diffuses at large time as ^(t) ^ 2Dt with a diffusion coefficient: 

an integral which converges at small s when the process is smooth since then c'(0) = 0. In fact, this 



formula, as well as (19 1 holds also for some processes with c"(0) = +oo, e.g. such as c(t) = e ^ at 



small T with a > 1, the main condition being that the small time singularity of C^(r) be integrable. 
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The convergence at large s should be discussed separately. Since D > J dsc'{s)'^ a necessary condition 
for convergence at large s is / dsc'{s)^ — J ^uj'^c{uj)'^ < +00. For a positive c(s) this is guaranteed by 
c(oo) — 0. For oscillating c(s), e.g. c(s) — cos{s)f{s) it requires |/(s)| to decay faster than l/y/s, or 



in Fourier, e.g. 
as: 



if c(w) ~ l/|w — then one must have 5 < 1. Interestingly, (20) can also be written 



D 



(21) 



where c(s) = sin6'(s) with 9{s) E [— 7r/2, 7r/2], ^(O) — 7r/2 and here 9{+oo) — 0. From there one 
sees that the strict criterion for convergence is that / ^uj'^cg{uj)^ < +00 where C6)(cj) is the Fourier 
transform of ce{T) — {9{0)9{t)), and is also positive (the variable 9 however is not Gaussian). 

Let us give examples of some of the non-generic situations where winding angle diffusion 
does not occur. The simplest case is the process — Cie'* ^ 
complex gaussian noises, which has correlation cu' — cos(i — t'). 



1^26 with ^1 and ^2 two i.i.d. 



One finds Cu(t) = — 



In sin 

2cos2 



and 



dr^ij) = 2r — tanrln(sin r), hence the winding angle grows faster than diffusive as ^{t) ^ t 
Consider next c(t) = Jo(t), i.e. c{uj) 
s and one finds superdiffusion ^(t) 



^(l — cj-^)(l — w^)^/^. The integral ( |20| is log-divergent at large 
^rlnr at large r. A range of superdiffusion can be obtained 
e.g. c(t) ^ T^''sm{T + tp) at large r yields ^{t) ^ t^"^^ for < 6 < 1/2. 

The above predictions are checked numerically in Fig. [3]in the time variable t, and as a parametric 
plot using c(t) in Fig. |4] for the diffusive and superdiffusive case. 




Figure 3. The angular velocity correlation function Cv{t) as a function of r for the c(r) used in 
Fig.l (the blue curve with the stronger decay), and for c(t) = Jo{t) (the second blue curve) together 
with the predictions of Eq. |18| (red curves). In the inset the winding angle variance ^(t), divided 
by r, is displayed for the same two cases. From the top, the first function (in blue) saturates to its 
diffusive value (red line), with D ~ 2.92 calculated from Eq. ( |20| . The second fu nct ion (in blue) is 
compared with the superdiffusion prediction "3?(t) = ^rlogr -|- 0.907t from Eq. |l9| . Both results 
are an average over 50 realizations 

Finally let us consider the stationary process Cair) = e'^~^/°^^^, which as a ^ converges to the 
non-smooth process C(r) = e"'^'. For any small a > the process is smooth and leads to diffusion. 
The diffusion coefficient however diverges as D « ^ ln(l/a) -I- ^ as a ^ 0. As we will see below the 
non-smooth limit process C(t) — e^'"^' is related to Brownian motion and leads to broad distributions 
of winding angle and to an infinite variance due to singular behaviour at short times. 
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c(r) c(t] 

Figure 4. Left panel: Parametric plot of the winding angle variance ^{t) (j/-axis) and c(r) (x-axis) 
for the c(r) used in Fig.l (blue curves) and the linear diffusion formula "I'(t) = 2Dt with D ~ 2.92 
as predicted from Eq. ( |20[ l (in red). The inset shows the angular velocity correlation function Cvij) 
as a function of c(r) for for the same choice of c(t) (blue curves) and the results of Eq. (red 
curves). Right panel: Parametric plot of the winding angle variance ^{t) as a function of c(t) 
for c(t) = Jo(r) (blue curves) and the asymptotic prediction 'I'(t) = ;^rlogr + 0.907r calculated 
from Eq. (19) (in red). In the inset the correlation function Ci,(r) is shown as a function of c(t) for 
the same c(r) (blue curves), together with the results of Eq. ( |18| (red curves). Both results arc an 
average over 50 realizations 



4-.2. non- stationary processes 

Wc now study non-stationary processes. Such process often occur in the context of aging or coarsening 
dynamics [H]. In fact they can, in some cases, be mapped onto a stationary process using the 
property of reparametrization of time: if the process Cu' has a winding angle (fit then the process 
Cg(()g(tM has a winding angle 4ig(t) for any positive monotonic function g{t). Note that Eq. (16 1 



and (17) have precisely the differential form required to satisfy this property. Hence for processes 
of the form cu' = c{g{t) — g(t')), the variance of the winding angle is immediately obtained as 

— ^{git) ~ 9{t')) where <&(t) is the variance for the stationary process c(r). Hence diffusion 
in $(t) ^ 2Dt implies ^tt' = 2(5(t) — g{t'))D at well separated times. One example, frequent in aging 
processes, is Cu' — f{t'/t) for t > t' . Then one can choose g{t) = Ini and c(s) = f{e~'^). To avoid 
divergences at short time difference one must have that /(x) = 1 — (1 — x)^ for x close to unity, with 
6 > 1. One example is f{x) = e~|i"^l . If this is the case, and provided convergence at large s, one 
finds that ^tt' ~ 2Dln{t/t') at well separated times, i.e. diffusion in the logarithm of time. Clearly 
the case 6=1 leads again to a non-smooth process and is discussed below. 

Among non-stationary processes, processes with stationary increments are of special importance. 
One such process is the so-called fractional Brownian motion (FBM), Cu' = \{t^'^ + (i')^'' ~ I* ~ ^T'*); 
with < /i < 1, which is the only member of this class which is also scale invariant. For ft, = 1/2 one 
recovers the standard Brownian motion. The FBM with h> 1/2 is sufficiently smooth for the above 
considerations to apply and one easily sees that the time change git) = In t and 

c(s) = cosh(/is) - 22''-i| sinh(s/2)|2'' (22) 

can be used, leading to diffusion for the winding angle in the variable g{t) = Int at large times, i.e. 
<^tt' - 2Dhlii{t/t') where = c'(s)V(l - c(s)2) diverges as /i 1/2+. 

The cases of the Brownian motion h = 1/2 and of the FBM for h < 1/2 require a separate 
discussion. Let us first recall what was found for the two dimensional BM with diffusion coefficient 
Dq i.e (ICtP) = 2Cft = 2DQt. At large t, for BM in the full plane the classical result [8J is that 
y — 2(j)t/\nt has a Cauchy distribution po(y) = -^{l + V^)^^ , hence the variance <^tt' cannot be defined. 
This broad distribution is regularized in presence of an absorbing (respectively reflecting) center in 
of radius R, where the distribution of y = 2(j)t / ln{DQt / R"^)) becomes PaIv) — 7r/(4cosh^(7ry/2) (resp. 
PB.{y) — l/(2cosh(7r?;/2))). In the two latter cases one has <i>tf/ ^ In^ (Dot / R^) at large t and fixed 
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t' (resp. <i>tt' ~ I In [DqI/ R'^)), however these are no more gaussian processes so the comparison with 
our results is not straightforward. We see below however that the process studied here yields a very 
similar result. Before we do so let us recall how the above scaling (j)t ~ Int for the (regularized) BM 
results can be understood from simple arguments. From the BM properties one easily obtains the 
stochastic equations for radius and angle (in Ito formulation) as drt — dBt + ^ and d^t — dBt/rt 
where, for Dq = 1, Bt, Bt are two independent unit Brownian motions (i.e. with (dBf) = (dBf) = dt). 
Diffusion in the winding angle is thus only possible if (l/r^) is bounded, as also nicely discussed in 
[T^ . If the Brownian can explore large distance then ^ ln{Dot/R^)/{Dot), where the small 

distance cutoff is also necessary, and one recovers the above non-diffusive behaviour (j)t ~ In t from the 
estimate = Dadt{l/rf) - dt\n{Dotl R^) /t. 



Can we make contact with our results, in particular can we also obtain from our formula ( 16 1 
the (regularized) BM scaling <^tt' ^ In^ t ? The answer is yes, but since we can only address smooth 
processes, we now consider the general smooth process with stationary increments: 

Cu' = \{f{t) + f{t')~f{t-t')) (23) 

with f"{t) = 2C2(t) in the notations of the Introduction, with /(O) = hence Cu = f{t)- The choice 
f(t)^t at large t corresponds to the random walk with a short time cutoff. Apart from the short times, 
it should look like the BM on large time scales. One example is = /J dt'rjt' where {rjtTjt') = ^e"'*"* ' 
then f{t) = i — 1 + e^* ~ /2 at short times. In general f{t) is an increasing function. Taking the 
large t limit at fixed t = t — t2 one finds: 

C.(t,i-r)«-^ln(/(r)//(t)) (24) 

This gives dt^tt' — 2 * drCvlt, t — t) k, ^y|^ /q°° dTf"{T) hence for the random walk /(<) ~ D^t 
at large t one finds: 

<^u' - lilntf (25) 

and one recovers the behaviour of the regularized BM. The prefactor, however, is different from both 
the absorbing and refiecting core, but is nicely anticipated from the simple argument presented above. 

The same calculation can be performed for the regularized FBM for h < 1/2, i.e. for the FBM 
random walk with f{t) ^ t^^ at large t with /i < 1/2 and f{t) smooth at small t. There one finds that 
the leading term above vanishes and one obtains: 

i.e. a much faster growth of the variance of the winding angle. 

Note that apart from the case of exact scale invariance f{t) = t^^ , there is no time 



reparametrization which allows to map the problem (23 1 with an arbitrary J{t) to a stationary process. 
And only for h > 1/2 the stationary process corresponding to the FBM, f{t) = t'^^, is smooth enough 



so that the present results can be used: note that in that case, the contribution (24) of the regime of 
fixed t — t2 is integrable and contributes only a constant to the winding variance, while the regime t2/t 
fixed (usually called aging regime when dealing with two time correlations) gives the main contribution, 
leading to the result ^tt' ^ 2Dfiln{t/t') found above. Conversely, the aging regime gives exactly zero 
contribution for the BM h — 1/2. Indeed then Cu' ~ t + f — \t ~ t'\, hence Cu' — /t for t > t' , and 
Cv{t,t') = for t > t' with a singularity at t = t' . Similarly, the aging regime is subdominant for the 
FBM random walk with /i < 1/2. 



4-.3. from Spitzer result to the winding of the stationary process c(t) = e 



-\r\ 



Let us return to the unrestricted BM motion for which Cu' = min(t,t') hence Ctt' = \/V jt for t > t' . 
Clearly because of the singularity sX t = t' in Cy(t,t') one cannot apply our formula for the winding 
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Figure 5. The v aria nce of the winding angle for t' = 1 as a function of t for the random walk, 
i.e. Cf^i of Eq. (|23| with f{t) = i — 1 + e~* (bottom blue curve) as compared to the asymptotic 
prediction of Eq.( |25[ l (corresponding red curve). The correlation functions for t' = 1 as a function 
of t for the FBM with h = 0.6 calculated using the equivalent stationary process \22\ with the time 
reparametrization s — > e* (top blue curve). The asymptotic diffusion prediction, = 2Di^—q^ logt 
where Eq. | |20| l gives -Dh=o.6 ~ 1-7 is also shown (top red curve). The results for the random walk 
required an average over ~ 10^ realizations. Inset:The rat io ^ ^^i / IrP t for the random walk, same 
data, showing the convergence towards 1/2 as predicted in ||25|l 



angle to this case. However the formula ( [T0| ) does hold for the BM and one has for integer n: 

^gm(0t = Fni^/tjt') (27) 

with the large t decay (at fixed f) (e™('^'*-"^'t')) _ ^'^+,^'' (^7^)"/^. 

Using the time change g(t) = In t one can now use the reverse correspondence and transport the 
Spitzer result [5] for planar BM to the stationary process c(r) = e~2l'^l which is not smooth and 
cannot be analyzed with the above methods. Hence the prediction for this process is that y — 2(I)t/t 
is distributed at large r with the Cauchy distribution pa{y). Its variance ^tt' is thus infinite at all 
times, as for the BM. 

An example of such a process is the Brownian motion or an ideal chain in an harmonic well, with 
C{u;) = l/(w^ + m^), hence in real time C(r) = j^^e"™'^' and c(r) = e""*!"^'. The distribution of 
the winding for such a confined Brownian is thus again the unit Lorentzian distribution for the scaled 
variable y = (pr /mr. An interesting generalization of a discrete version of this model to a chain, was 
studied in Ref. [17) . with the result that, again, each monomer sees a Lorentzian winding. 

It is interesting to now consider a smoother variant of this model, i.e. an ideal chain with a small 
curvature energy in a harmonic well, described by C{lo) — (w^ + m^)^^ — {uP' + A'P)^^ k, 1/(0;^ /M^ + 
+ m?) at large M ^ m. The decay in the time domain, c(t) — (Me^™^ — me^^^'^)/{M — m) is 
now smooth at small times, and one finds that the winding angle recovers now a finite variance and is 
diffusive <^tt' ^ 2Dt with a diffusion coefficient D - iMln(4.2/M). 
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5. Algebraic area enclosed 

Finally we can study the algebraic area At enclosed by the process, which satisfies At = \{£,t S,t ~ ^t) ■ 
Its variance is extracted from Go^Q^t.t' above and one easily finds that 

CA{t, t') = {AtAt.) = lidt'/^Ct^t' - 4' "^4' (28) 

For a smooth stationary process one finds Ca{t) — ~^C" {t)C{t) + ^C'{t)'^, and dr{A'^{T)) = 
2 loCAir) = 2 /J'C7'(s)2ds - C"(r)C(T) and one finds the diffusion result {A^{t)) ~ 2Dat with 
Da — /q°° C'{s)^ds. Let us consider now the above process with stationary increments. Note that the 
time reparametrization is useless here. For the random walk f{t) ~ D^t one finds {A^) ^ Dot^/A at 
large t. This is larger than the result for Brownian paths constrained to come back to their starting 
points (loops) obtained in Ref. [H]. This is well confirmed by our numerics displayed in Fig. [6] where 
the result for a stationary process, which instead exhibits only diffusive growth of the area, is also 
shown for comparison. 




Figure 6. The variance of the algebraic area 'I'^ = {[At — ^t']^) with t' = 1 as defined above 
Eq.| |28| for : (i) the random walk, i.e. Ctti as in Fig. 5 (blue curve on top at large time) and the 
asymptotic prediction ~ (corresponding red curve) (ii) the result for the stationary process c(t) 
of Fig.l as a function of t = t (bottom curve) and the asymptotic prediction 2D At with Da = 37r/8 
(corresponding red curve). In the inset the ratio ^a/^^ is plotted for the random walk as a function 
of t, and shows saturation towards the predicted prefactor 1/4. 



6. Conclusion 

We have computed here the angular velocity correlation of a very general smooth gaussian process 
in the plane. This allowed us to obtain a simple closed formula for the diffusion coefficient of the 
winding angle valid for most such stationary processes. Our formula also extends to non-stationary 
processes, and has allowed us to obtain the three main behaviours (i) diffusion in the logarithm of 
time for sufficiently smooth fractional Brownian motion (ii) the square of logarithm in time (25) for 
the winding for the usual random walks related to the Brownian motion (iii) power law growth in 
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time ( 26 ) for the winding angle of the random walks which provide a regularization of the non-smooth 
fractional Brownian motion. 

Various extensions of the present calculations are left for the future. These include: higher 
moments and distributions of area and winding, winding around a point different from the origin, or 
around several points, winding conditioned to closed paths, most general 2D gaussian process including 
non zero average, and finally, devising methods to account perturbatively for non gaussian effects. It 
would also be interesting to study persistence effects such as the probability that the winding angle 
never crosses zero, or to compute the winding for loops, i.e. conditioning the process to return to its 
starting point. 
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